Multiple ordered phases in the filled skutterudite compound PrOs4Asi2 
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Magnetization, specific heat, and electrical resistivity measurements were made on single crystals 
of the filled skutterudite compound PrOs4Asi2. Specific heat measurements indicate an electronic 
specific heat coefficient 7 ~ 50—200 mj/mol at temperatures 10 K < T < 18 K, and ~ 1 J/mol 
for T < 1.6 K. Magnetization, specific heat, and electrical resistivity measurements reveal the 
presence of two, or possibly three, ordered phases at temperatures below ~ 2.3 K and in fields below 
~ 3 T. The low temperature phase displays antiferromagnetic characteristics, while the nature of 
the ordering in the other phase(s) has yet to be determined. 

PACS numbers: 71.27.+a, 75.30.Kz 



The filled skutterudite compounds, with the formula 
MT4X12 {M = alkali metal, alkaline-earth, lanthanide, 
or actinide; T = Fe, Ru, or Os; and X = P, As, or 
Sb), display a wide variety of strongly correlated elec- 
tron phenomenaiSi^ Of particular interest are the Pr- 
based filled skutterudites. The physical properties of 
these compounds are dominated by the ground state and 
low lying excited state of the Pr'^^ ion in the crystalline 
electric field (CEF), and the hybridization of the Pr 4f- 
states with the ligand-states of the surrounding Sb ions 
that compose the atomic cage within which each Pr"^+ 
ion resides. A variety of correlated electron phenomena 
have been observed in the Pr-based filled skutterudites: 
conventional (BCS) and unconventional superconductiv- 
ity, magnetic order, quadrupolar order, metal-insulator 
transitions, Kondo phenomena, heavy fermion behavior, 
and non-Fermi liquid behavior. A prime example of this 
diversity is seen in the compound PrFe4Pi2, which un- 
dergoes a transition to an ordered state below 6.5 K that 
was originally thought to be antiferromagnetic (AFM) in 
nature,"* but was later identified with antiferroquadrupo- 
lar (AFQ) order The suppression of this AFQ ordering 
in magnetic fields results in the formation of a heavy 
Fermi liquid state above the quadrupolar quantum criti- 
cal point (QCP)&2i& Below 1 K, a possible high field or- 
dered phase between 8 T and 12 T has also been observed 
in a limited angular range around the [111] directioniS 
One of the most interesting Pr-based materials is the 
compound PrOs4Sbi2, which exhibits unconventional su- 
perconductivity, involving heavy fermion quasiparticles, 
below a critical temperature Tc = 1.85 KjiSiiiii^ The 
superconducting state breaks time reversal symmetry^ 
appears to consist of several distinct superconducting 
phases fi2iiiii& and may have point nodes in the energy 
gg^p^i4ji6ji7 ■^]^[Iq between 4.5 T and 16 T and below 



~ 1 K, an ordered phase is observed that has been iden- 
tified with AFQ orderiSiiLiSiia This suggests that the 
extraordinary normal and superconducting properties of 
PrOs4Sbi2 may be associated with the proximity to a 
field-induced quadrupolar QCP, in accord with previous 
conj ectures j -^°i**i*^ 

Unlike the Pr-based filled skutterudite phosphides and 
antimonides, the arsenides have not been investigated in 
much detail. These materials provide an opportunity 
to discover other examples of strongly correlated elec- 
tron behavior displayed by the Pr-based filled skutteru- 
dites. In this report, we present an investigation of the 
physical properties of PrOs4Asi2 single crystals. These 
measurements reveal that PrOs4Asi2 is a Kondo lattice 
compound with an enhanced electronic specific heat co- 
efficient 7 of ~ 50 — 200 mJ/mol at temperatures 
10 K < T < 18 K, and ~ 1 J/mol for T < 1.6 K. 
Multiple ordered phases are observed below 2.3 K and 
3.2 T, at least one of which appears to be AFM in na- 
ture. 

Single crystals of PrOs4Asi2 were grown from elements 
with purities > 99.9% by a molten metal flux method at 
high temperatures and pressures, the details of which will 
be reported elsewherejSS. After removing the majority of 
the flux by distillation, PrOs4Asi2 single crystals of an 
isometric form with dimensions up to ~ 0.7 mm were col- 
lected and cleaned further in acid to remove any Pr impu- 
rity phases from the surfaces. X-ray powder diffraction 
measurements were performed with a Rigaku D /MAX B 
X-ray machine on a powder prepared by grinding several 
single crystals along with a high purity Si(8N) standard. 
The crystal structure of PrOs4Asi2 was determined by 
X-ray diffraction (XRD) on a crystal with a regular oc- 
tahedral shape and dimensions of 0.12 x 0.12 x 0.12 mm. 
A total of 5158 reflections (583 unique, Rint = 0.1249) 
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TABLE L Atomic coordinates, displacement parameters, and occupancy factors for PrOs4Asi2. \Jeq is defined as one third of 
tlie trace of tlie ortliogonalized Vij tensor. 
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FIG. 1: (a) Magnetic susceptibility Xdc vs T, indicating a 
magnetic transition at 2.3 K. (b) Inverse magnetic suscepti- 
bilility Xdc '^^ ^- The line represents a Curie- Weiss law, fit 
to data above 50 K. The inset shows the low-T behavior of 
Xdci"^) s^long with a different Curie- Weiss law fit (see text). 



were recorded and the structure was resolved by the fuU 
matrix least squares method using the SHELX-97 pro- 
gram with a final discrepancy factor Rl = 0.0462 (for 
I > 2ct(I), wR2 = 0.1124),aii22i 

Magnetization M measurements for temperatures 
0.4 K < T < 10 K and magnetic fields -5.5 T < H < 
5.5 T were made in a ^He Faraday magnetometer. A 
mosaic of 27 crystals, with a total mass of 27 mg, was 
measured with the [111] axis of the crystals nominally 
aligned parallel to H. Measurements of M{T,H) were 
also performed using a Quantum Design Magnetic Prop- 
erties Measurement System (MPMS) in H up to 5.5 T 
for 1.7 K < T < 300 K. Four- wire ac electrical resistivity 
p{T) measurements were made for 2 K < T < 300 K in a 
Quantum Design Physical Properties Measurement Sys- 
tem (PPMS) using a constant current of 1 — 10 mA, and 
for 0.06 K < T < 2.6 K in a ^He-^'He dilution refrigerator 
with a constant current of 300 /iA. Measurements of p{T) 
at various pressures P up to ~ 23 kbar were made down 
to 1 K in a ^He cryostat using a BeCu piston-cylinder 
clamp. Specific heat C(T) measurements on PrOs4Asi2 
were performed for 0.6 K < T < 35 K at UCSD in a 
semi-adiabatic '^He calorimeter on a collection of single 
crystals with a total mass of 49 mg. Additional C(T) 



FIG. 2: Magnetization Af vs H isotherms at several values of 
T. Features in M{H) at constant T, with Hi and H2 defined 
by a kink in M{H), are illustrated for T = 0.425 K. The inset 
shows M{H) at higher temperatures. 



measurements for 0.4 K < T < 20 K were performed at 
Lawrence Livermore National Laboratory (LLNL) in a 
Quantum Design ^He PPMS. 

Analysis of the X-ray powder diffraction pattern in- 
dicate single phase PrOs4Asi2 with no major impu- 
rity peaks. A unit cell parameter a — 8.5319(11) A 
was determined, in excellent agreement with an earlier 
measurement^^ of a = 8.5311(3) A. Single crystal struc- 
tural refinement shows that the unit cell of PrOs4Asi2 
has the LaFe4Pi2-type structure (Im3 space group) with 
two formula units per unit cell, and a = 8.520(1) A. 
Other crystal structure parameters are summarized in 
Table The displacement parameter Ue^ represents the 
average displacement of atoms vibrating around lattice 
positions and equals mean-square displacements along 
the cartesian axes. The displacement parameters deter- 
mined for PrOs4Asi2 show behavior which is typical of 
the lanthanide filled skutteruditesi^^iS^ Table also indi- 
cates that the Pr site in PrOs4Asi2 can be assumed to 
be fully occupied since there is a ~ 2% uncertainty in the 
occupancy factors. 

The dc magnetic susceptibility Xdc(T), measured in 
H = 0.05 T along the [111] direction, for 0.4 K < 
T < 10 K, is shown in Fig. Q^a). The clear peak in 
the Xdc{T) data suggests an AFM transition at T = 
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2.3 K. Figure Ulb) displays xlc vs T in F = 0.1 T be- 
tween 2 K and 300 K. Above 50 K, the data are well- 
described by a Curie- Weiss law with an effective mo- 
ment /icff = 3.81/^B/f-u. and a Curie- Weiss temperature 
Bcw = -25.8 K. Below 20 K, x^^^T) can be described 
by a different Curie- Weiss law, with /icff = 2.77/iB/f-u. 
and 0CW — 0.44 K [Fig. ^b) inset]; this behavior is 
consistent with crystalline electric field (CEF) effects, 
as will be discussed later. The value of ^cS derived 
from the high temperature (above 50 K) Curie- Weiss 
fit is larger than the theoretical Pr'^^ free-ion value of 
HeS = 3.58/iB/f.u. A similar /ies = 3.84//B/f-u. was 
also found for PrRu4Pi2.^^ In contrast, PrOs4Sbi2 dis- 
plays an effective moment of 2.97/iB/f-u.fii while both 
PrOs4Pi2 and PrFe4Pi2 display effective moments close 
to the Pr''+ free-ion value. '''^^ The larger-than-expected 
moment may be due to an unaccounted-for constant 
Pauli paramagnetic contribution. When a constant con- 
tribution xo of 5.5 X 10^"* cni"^/mol is assumed, which 
is consistent with a value of 5.16 x 10^^ cm'^/mol found 
for LaFe4Pi2ri the effective moment estimated from a 
Curie- Weiss law fit approaches the Pr^+ free-ion value 
with Qcw = —12.5 K. The values of ^cS and Gcw for 
the low temperature (below 20 K) Curie- Weiss law fit re- 
main unchanged when the Xdc{T) data are corrected for 
this value of xo • 

Low-T M{H) data taken with H parallel to [111] are 
shown in Fig.|21 An infiection point at low H is evident, 
along with two kinks in M{H) (denoted Hi and H2) 
which are suppressed toward lower field with increasing 
T. The inset of Fig. |21shows the evolution of M{H) from 
higher temperatures down to the ordering transition tem- 
perature. No magnetic hysteresis was observed. Features 
corresponding to Hi and H2 are found in M{T) in con- 
stant H. With increasing 77, a peak in M{T) [labeled as 
Ti in Fig. ^a)] shifts towards T = K, while a kink in 
M{T) at higher T (denoted T2) becomes apparent, splits 
off from Ti, and eventually also moves towards T = K. 

Specific heat C(T) measurements on PrOs4Asi2 in zero 
field (Fig. |3Jl display two clear features: a slight shoul- 
der at T2 ~ 2.3 K, on the high temperature side of a 
well defined peak at Ti 2.2 K. Zero-field measure- 
ments made at LLNL show similar behavior [Fig. 13 inset 
(b)], except that in the LLNL experiment, the peak at 
Ti was further resolved into two distinct features sep- 
arated by ^ 0.02 K. This additional feature may in- 
dicate an additional transition or it may be due to a 
variation in composition among the measured crystals. 
At the lowest temperatures, there is a slight upturn in 
C{T) associated with a nuclear Schottky anomaly. The 
low-T specific heat can be described by four contribu- 
tions, C{T) = Cn{T) + Co{T) + Ci{T) + C^{T), which 
are the nuclear Schottky, electronic, lattice, and mag- 
netic terms, respectively. At higher temperatures, Cm{T) 
and Cn{T) w 0, and the specific heat is dominated by 
the electronic and lattice terms. In Fig. 0{a), a fit of 
C/T = 7 -t- /3T2 in the range 10 K < T < 18 K re- 
veals an enhanced 7 « 200 mJ/mol and, from (3, a 
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FIG. 3: Specific heat C divided by T vs T (UCSD data). 
Inset (a) shows the transitions at Ti and T2 in more detail 
(UCSD data). Inset (b) shows the apparent resolution of the 
transition at Ti into two transitions at Ti and Ti' (LLNL 
data). 
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FIG. 4: (a) Zero field C/T vs T^, with a linear best fit yield- 
ing the enhanced 7 and value of 0_d given in the figure, (b) 
Magnetic entropy Sm determined from the specific heat af- 
ter the zero field lattice and electronic contributions to the 
specific heat were subtracted. 



Debye temperature Qd ~ 260 K. To analyze the spe- 
cific heat in the ordered state, the three terms of the 
specific heat after subtraction of the lattice contribution 
(using w 260 K), AC{T) = Cn + Ce + Cm, were 
assumed to have the following forms: Cn(T) — A/T"^, 
Cc{T) = 7T, and Cm(T) = BT". Using these equations, 
fits of ACT2 = A + 7T3 + BT"+2 were performed for 
T < 0.6ri. The best fit results in a nuclear Schottky 
coefficient A = 128 mJ K/mol and an exponent in the 
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FIG. 5: Electrical resistivity p vs T from 2 K to 300 K. The 
inset shows p{T) from 0.06 K to 25 K. 



power-law term of n = 3.2 ± 0.1, close to the expected 
value of n = 3 for AFM spin waves. The electronic 
specific heat coefficient 7 is strongly enhanced, with a 
value near 1 J/mol K^. A temperature dependent 7 is 
not uncommon and has been observed in numerous com- 
pounds that exhibit heavy fermion behavior. ■^^ The value 
of 260 K for is consistent with other filled skutteru- 
dite arsenides with typical values ranging from 230 K 
(LaRu4Asi2) to 340 K (PrRu4Asi2).22i2i The magnetic 
entropy SmiT) was calculated from the specific heat by 
integrating Cm (T) /T with temperature [Fig. Elb)] , re- 
sulting in an entropy release of 90% i? In 3 that levels off 
by 8 K, a value consistent with a triply degenerate ground 
state. The reduction of the entropy from a full i? In 3 is 
presumably due to hybridization of the localized 4f and 
itinerant electron states, as is observed for PrOs4Sbi2, 
where entropy is apparently transferred from the local- 
ized 4f states to the conduction electron states. 

The electrical resistivity p(T) was measured in the 
range 0.06 K < T < 300 K (Fig. EJ. Kondo lattice be- 
havior is evident, as p(T) decreases with decreasing T 
down to Tniin ~ 16.5 K and then increases until 5.5 K, 
where the resistivity drops due to the onset of the or- 
dered state. Measurements in high magnetic fields can 
be analyzed in terms of a single ion Kondo model for an 
impurity effective spin of 1, yielding an estimate of Tk 
on the order of 1 K«S2, A small Tr- with respect to the 
ordering temperature(~ 2.3 K) prohibits the observation 
of dependence of p{T) . Measurements of p{T) under 
P up to 23 kbar between 1 K and 300 K revealed no 
significant P dependence. 

As mentioned earlier, features in the XdciT) data re- 
flect CEF effects. Although the true crystal symmetry 
of the filled skutterudites is tetrahedral Th, it is a close 
approximation in zero field to consider cubic Oh sym- 
metry. In a CEF with Oh symmetry, the ninefold de- 



FIG. 6: (a) Schematic representation of the PrOs4Asi2 mag- 
netic field - temperature {H — T) phase diagram. The vertical 
and horizontal lines represent isothermal M{H) and constant 
field M{T) measurements as depicted in (b) and (c). (b) Mag- 
netization M vs H at 0.425 K, with features at H\ and H2 
used to define the phase diagram, (c) Magnetization M vs T 
at 1.5 T, along with two clear features at T\ and T2, used to 
contrast the phase diagram, (d) The derivative dM/dH as a 
function of H indicates two features used to define H\ and H2 
in Fig. inib) . The dotted lines represent estimated transition 
widths for H\ and i?2. 



generate Pr'^^ J = 4 multiplet is split into a Fi singlet, 
a Pa doublet, and r4 and F5 triplets. The low-T cal- 
culated effective moment of 2.77/iB/f.u. is indicative of a 
Fs ground state, which has a theoretical effective moment 
of 2.83^B/f-u. In addition, M = 1.68/iB/f-u. at 2 K and 
5.5 T, a value much less than the full 3.2/iB/f-u. saturated 
moment expected for a Pr'^^" free ion, but more consis- 
tent with the 2.0/iB/f.u. expected for a F5 ground state. 
Reasonable CEF fits to the Xdc {T) — xo data (not shown) 
can be made using a CEF splitting as low as ~ 15 K be- 
tween a Ps ground state and a Fi first excited state. If 
there were a large Schottky contribution due to this CEF 
splitting in the temperature range of the high tempera- 
ture fit to C(T), the values of 7 and Qd determined from 
the fit would be altered. Based on the CEF fits to the 
Xdc{T) — xo data, the worst case scenario (which would 
lead to the largest possible change in 7) occurs for a 
splitting of 23 K. The resulting Schottky contribution to 
the specific heat would result in a reduced 7 (50 mJ/mol 
K2) and a shghtly lower Gd (246 K). Since the CEF 
contribution to C(T) is not significant below 1.6 K for 
splittings above 15 K, the large value of 7 (~ 1 J/mol 
K^) determined from the low temperature fit to C{T) is 
not affected by potential CEF effects. The various ac- 
ceptable CEF splitting schemes do not significantly alter 
the magnetic entropy release of 90% i? In 3 by 8 K, which 
remains in support of a triplet ground state. 

A schematic representation of the {H — T) phase dia- 



5 



4 



2.5 





PrOs As 

4 12 


o 


J 




Magnetic Phase Diagram 


• 

□ 


1 

\ 

Hi 
H 

2 


- 




■ 




1 — Of 

1 — ax 















0.5 1 1.5 2 2.5 



T(K) 



FIG. 7: Magnetic field-temperature {H — T) phase diagram 
for PrOs4Asi2, as defined by features in the magnetization 
(Fig. HKcircles and squares). The hues are guides to the eye. 

gram for PrOs4Asi2 (Fig. |7I), is shown in Fig. [Sfa). The 
vertical line labeled "b" in the schematic phase diagram 
represents an M{H) measurement at T = 0.425 K shown 
in Fig. Elb). The points Hi and H2 along line b are 
features in the M{H) data, which are defined using the 
derivative dM/dH shown in Fig. Eld) as a function of 
H. The horizontal line labeled "c" represents an M{T) 
measurement at H = 1.5 T (Fig. (Sfc)). The features at 
Ti and T2 along c are defined from the M{T) data shown 
in Fig. ini^c). Similar features in the other M{H,T) mea- 
surements were used to generate the full H ~T phase di- 
agram for PrOs4 Asi2 shown in Fig. |7| The error bars are 
rough estimates of the transition widths, based upon the 
widths of features in the derivatives of M{H) (Fig.E^d)) 
and M{T) (not shown). It may be inferred from the fea- 



tures in the UCSD C(T) data that there are two phase 
transitions. Thus, PrOs4Asi2 appears to have an AFM 
ground state, a second ordered state at intermediate T 
and H ^ and is paramagnetic at high T and H . 

The AFM nature of the ground state is supported by 
several observations. The peak in Xdc{T) is a traditional 
indication of the onset of AFM, which is also suggested 
by the lack of hysteresis in M{H), and is compatible with 
a CEF-split T5 magnetic triplet ground state. Further- 
more, the suppression of the ground state order with H 
is consistent with AFM order. The nature of the second 
ordered phase is difficult to establish from the results of 
this study, but in light of the low-T behavior of other Pr- 
based filled skutterudites such as PrFe4Pi2ii^ magnetic 
or quadrupolar ordering seem likely possibilities. Recent 
neutron diffraction experiments on PrOs4Asi2 confirm 
that the low field ordered phase is antiferromagnetic.^^ 

To summarize, magnetization and specific heat mea- 
surements on PrOs4Asi2 indicate the existence of at 
least two distinct phase transitions at temperatures be- 
low 2.3 K and in fields below 3 T. The magnetic behavior 
of the low temperature phase is consistent with AFM or- 
der. The electrical resistivity displays a sharp decrease 
due to the onset of ordering and has features consistent 
with Kondo lattice behavior. In the paramagnetic state, 
specific heat measurements indicate an electronic specific 
heat coefficient with 7 « 50 — 200 mJ/mol K^. In the 
antiferromagnetic state, analysis of C(T) data yields a 
very large 7^1 J/mol at low T (T < 1.6 K). 
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DE-FG52-03NA00068. 
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